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Abstract 

We make an adaptation of laser modelling equations to describe the behavior of a 
phonon laser (saser). Our saser consists of an AlGaAs/GaAs double barrier het- 
erostructure designed to generate an intense beam of transversal acoustic (TA) 
phonons. To study our system, we begin with a Hamiltonian that describes the de- 
cay of primary longitudinal optical phonons (LOi) into secondary {LO2) and TA 
(LOi LO2 + TA) and its inverse process (recombination). Using this Hamilto- 
nian, a set of coupled equations of motion for the phonons is obtained. We also 
consider the interaction between the phonons and its reservoirs. These interactions 
are introduced in the equations of motion leading to a set of coupled Langevin 
equations. In order to obtain an expression to describe our saser we apply, in the 
Langevin equations, an adiabatic elimination of some variables of the subsystem. 
Following the method above we obtain the value of the injection threshold for the 
operation of our phonon laser. At this threshold occurs a phase transition from a 
disordered to a coherent state. It is shown that it is not necessary a big "optical" 
pumping to get a sasing region. 
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1 Introduction 

Sasers are devices capable of producing beams of coherent phonons. Besides 
the system considered here, several other kinds of such devices were proposed 
and constructed [1-7]. 

In recent works [8-13] the kinetics and the dynamical properties of the double 
barrier phonon laser were discussed. In this paper we study the system from 
a point of view close to that used in quantum optics to study the onset of 
the lasing regime. In spite of the fact that the system is far from equilibrium 
and thus concepts like temperature and chemical potential can not be applied, 
the process is analogous to a phase transition [14,15]. This transition is not 
governed by temperature but an energy flux (in this case by the injection 
rate of primary longitudinal optical LOi phonons). We were able to obtain 
the potential that describes the states of our system and consequently, the 
threshold value for the coherent transversal acoustic {TA) phonon emission. 

The paper is organized as follows. A description of the saser device is presented 
in section 2. The multimode case is analyzed in section 3. In section 4 the 
potential that describes our system is obtained. Finally, in section 5 we present 
our results and conclusions. 



2 The system 

The phonon laser studied here consists of an AlGaAs/GaAs double-barrier 
heterostructure (DBH). This device was designed in such a way that, for a 
small applied bias V, the difference Ae between the first excited level £i and 
the ground state £o in the well is less than the longitudinal optical (LOi) 
phonon energy huoi. For a greater bias the resonant condition Ae ~ huji is 
achieved and the electrons begin to decay to the ground state by emitting 
primary LOi phonons. For an Al concentration greater than 0.25 [16] or 0.3 
[17], these phonons are confined inside the well (they can be also absorbed by 
exciting electrons from Eq to £i). 

The process described above, acts in parallel with the decay of primary LOi 
phonons due to anharmonicity. One of the products of this decay is a secondary 
longitudinal optical phonon {LO2), the other is a TA phonon [18]. The LO2- 
TA pair is produced by stimulated emission. Therefore these TA phonons 
could be coherent and form a beam that it is called saser by analogy with a 
laser. 

The DBH saser beam could be applied to acoustic nanoscopy, phonoelectronics 
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and phonolitography. This was discussed in more detail in a previous article 
[8]. 

The potential profile and the level positions at the resonant condition are 
shown in figure 1. 
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Fig. 1. Potential profile and energy levels at resonant condition, bi, br, d and ep are 
respectively the barriers and well widths and the Fermi level at the emitter. 



3 The multimode case 



The fundamental mechanism in this saser is the generation of the TA phonons 
due to the decay of primary LOi phonons. Therefore, the other processes (i.e., 
emission-absorption of the LOi phonon through electron transitions and the 
LO2 phonon decay) are secondary and are taken into account just as a result 
of the interaction of each phonon with its reservoir. 

As the TA phonon emission results from the decay of LOi phonons, they could 
be emitted in many modes, corresponding to the different values of its wave 
vector q3 parallel to the interface. 

The number N^^^ of wave vectors parallel to the interface can be estimated as 
A^qg ~ So/a^, where is the device area and a is the lattice parameter. For 
the values used here. So — 0.5 • 10~^ mm^ and a — 2.825 A, N^^ is of the order 
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of 7 • 10^ The total number of emitted TA phonons is ~ 2 ■ 10^ [13], which 
means that, at the beginning of the phonon emission, the number of phonons 
per wave vector qs is ~ 3-10^®. As the apphed bias increases, the emission in a 
particular mode called qg grow slowly until a certain threshold value V — Vth 
is reached. For values of V greater than Vth, the phonons in that mode are 
emitted in a coherent way implying that the population nqO is the only one 
macroscopically non zero. In this section, the possibility of macroscopically 
coexistence of more than one mode is studied. 

3.1 The Hamiltonian 

The Hamiltonian that describes our system can be written as: 

qi 12 qs 

^7 E Kb^K+Kblbl), (1) 
qiiq2,q3 

where qi, q2 and qs are respectively the LOi, LO2 and TA wave vectors 
parallel to the interfaces, and must satisfy qi = q2 + qs, {i = 1,2,3) are 
the energies of each kind of phonons and 7 is the phonon-phonon interaction 
coefficient. 

3.1.1 The Langevin equations 

The equations of motion for the operators riqj (riq. = 6^.6q.), b^^, b^^ are ob- 
tained from the Hamiltonian (1). We remark that the saser system is coupled 
to reservoirs. The LOi phonons are coupled through the emission and absorp- 
tion processes to the electron system. In the case of the LO2 phonons, they 
are coupled to the other phonon modes to which they decay, whereas the TA 
phonons are connected to the other TA phonon modes outside the well. Thus 
the total Hamiltonian includes, besides the terms shown in equation (1), the 
interaction with the reservoirs. By using the methods described in references 
[19,20], the reservoir coordinates can be eliminated leading to the following 
Langevin equations 

n E {blKK - Kblbl} + Gqi - rmq, + A, (2) 

q2,q3 

(i^2 + K2)&q2 - «7 E ^qi^is + -^2, (3) 
qi-qs 



dn, 



qi 



dt 

db, 



q2 



dt 
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^ = -{tUJs + K^,)b^, - «7 E KA^ + ^3, (4) 

qi,q2 

where Gq^ and Fi are the emission and absorption rates for the LOi phonons 
due the interaction with its reservoir (i.e., the electrons), K2 is the decay 
rate of LO2 phonons, k^^ is the escape rate of the TA phonons and Fi {i = 
1,2,3) are the fluctuation forces. The TA phonons escape rate is different for 
each mode, thus depending on q^. As we are not interested in studying the 
statistical properties of the saser, in the following the fluctuation forces will 
not be considered. 

In order to obtain only one equation for the TA phonons, we make the adi- 
abatic approximation in equations (2) and (3). First we replace 6qj (t) = 
^qi (t) and {t) = feqg {t) in (3). Secondly, we integrate (3) ex- 

phcitly to get 

t 

K = E / e-[^(— --3)+«.](t-r) {hJl^^dr. (5) 



Assuming that the relaxation time of the phonons is much smaller than 
the relaxation of 6q^ and 6^^, we can take ^qi&qg out of the integral in equation 
(5) obtaining 

K = -^7 E (6) 

qi,q3 -^^ + >^2 



where u — uji — UJ2 — oj^- 

In a similar way, we obtain for nq^ 

^q.=^(l-^E^q3)> (7) 



where 



Finally, we replace (6) and (7) in (4) and we get, after some manipulations, 
an equation for the TA phonon populations 

djTi 

^ = (Qni-2«q3)nq3, (8) 
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where rii — Y, ^qi ■ 
qi 

Considering as an example only two modes, and ql, we get for the stationary 
case of equation (8): 



If both modes would be present, (ng 7^ and nl ^ 0) both parenthesis in 
(9) and (10) would be zero. This is a contradiction (because k'^ 7^ kI). This 
contradiction only can be solved if just one mode is present and the other 
one has died out. Therefore, only a single mode, the one with largest lifetime 
inside the well and closest to resonance, survives. The former analysis also can 
be done quite rigorously for many modes. 

When the LOi decay, the TA phonons are emitted in the [111] direction [21]. 
The device is grown in this direction, such that the TA phonon beam will make 
multiple reflections between the walls of the well, implying in large lifetime for 
these phonons. Indeed, for a small applied bias, the TA phonons are emitted 
with any wave vector qg. In this regime, the device works as a sound emitting 
diode (SED). With a further increase in the applied potential, when a certain 
value of V is reached, the TA phonon distribution on qg become more and 
more sharp giving place to a phase transition. In this case grows suddenly 
from about 10~^ to about 10^. Thus, a great number of TA phonons would 
have only one mode: qs = (which imply qi = q2). This mode will survive 
and slave the others, making the system to self-organize and emit just in this 
mode. 



4 The single mode case 

In order to get the threshold for which begins to be macroscopically non 
zero, we must first obtain an equation for the amplitude of the single mode 
63 that slaves the system. To do that, we make the same approximations as 
described in the previous section. The equation obtained for 63 is 




(9) 
(10) 




(11) 



where d = E G, 
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The right-hand side of (11) can be obtained from the potential: 



Uc {ht) 



a 
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;i2) 



with 



a 



1^3 - 



r 



ri/«2 



r 



ri/«2. 



such that 



d6| 
dt 



dbl 



(13) 



The expressions for the parameters Gi and Fi are [13] 

Ti^w (No - Ni) , Gi = wNi 



where w is the electron transition rate and Nq and A^i are respectively the 
electron populations of the fundamental and excited levels. 

The shape of the potential (12) is similar to that obtained in the Ginzburg- 
Landau phase transition theory. For a > the potential has a minimum at 
^3 = (^3 = l&sl^)- In the case of a < 0, the minimum is obtained for 
^3 7^ 0. The transition from « > to a < is due to the competition between 

2 

processes of gain, with a rate Gi^^ and losses, with a rate k^. The threshold 
is attained when a = 0, i.e., Gi = ^^^^^ ■ 



5 Results and Conclusions 



The parameters used in our calculation are: the escape rate of TA phonons 
^3 = 0.05 ps~^, estimated from its group velocity near the L point of the 
Brillouin zone, the decay rate of the LO2 phonons K2 — 20 ps"^ (by assum- 
ing a small stimulated decay) and the phonon-phonon interaction coefficient 
7 = 0.5 ps~^ calculated in reference [22] following the work of Klemcns [23]. 
The decay rate 70 obtained from this value is in agreement with the exper- 
imental result of Vallee [18]. The values of Gi and Fi depend on the LOi 
phonon emission rate w and on the electron populations A^o and A^i that were 
calculated in reference [13] by solving the kinetic equations. 

The total generation rate of the LOi phonons Gi, depends on the applied 
bias V (i.e., it turns out from electron transitions that are V dependent [13]). 
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Fig. 2. The potential Uc as a function of the TA phonon amplitude for different 
values of the applied bias. 

In figure 2, the potential (12) is plotted as a function of the TA phonon 
amplitude for different values of the applied bias V. As it can be seen, for 
V > Vth = 59.4 mV this potential has different minima for different values of 
V. These minima represent the stable points 63 of the system. 

In figure 3 the stable points are plotted as function of V. The sohd hne rep- 
resents the n| calculated in this work. The dashed line is the TA phonon 
population calculated in a previous work from the kinetic equations [13]. 

Figure 4 shows a comparison in the region where the TA generation begins 
(dashed region in figure 3). From this comparison two conclusions arise. First, 
as it can be seen (solid line), the TA coherent emission does not start until 
a certain threshold of the applied bias is reached. At Vth the LOi emission 
rate Gi is of the order of 279 ps~^ and the total TA phonon population 
is ~ 4658, whereas — 0. The TA phonons generated for a bias greater 
than Vth s-re coherent. Secondly, our previous calculations [13], even without 
considering the self-organization process, gave results that agree quite well 
with those presented here. This results show that the coherent emission occurs 
immediately after the condition of LOi phonon emission is reached. That 
means that it is not necessary a big "optical" pumping to get a sasing region. In 
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Fig. 3. The solid line shows the result for the coherent mode population obtained 

here (the minima of the previous figure), whereas the dashed line shows the total 
TA phonon population obtained from the kinetic equations [13]. 

summary, we show in this work that the TA phonons, produced by stimulated 
emission, constitute a coherent phonon beam in the parameter region discussed 
above. 
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Fig. 4. The dashed region in figure 3 is shown in detail. We can see that the TA 
emission process begins incoherent. It becomes coherent when the injection rate of 
phonons LOi is Gi > 
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